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Abstract 

Low frequency longitudinal wave instability in magnetized high correlation 
dusty plasmas is investigated. The dust charging relaxation is taken into 
account. It is found that the instablity of wave is determined significantly 
by the frequency of wave, the dust charging relaxation, the shear viscosity 
and viscoelastic relaxation time, the coupling parameter of high correlation 
of dust as well the strength of magnetic field. 

PACS number(s): 52.25.Zb, 52.25.Vy, 52.35. Mw 

Low frequency longitudinal volume dust waves have been extensively studied for weakly 
coupled and unmagnetized dusty plasmas [1—5]. Recently the corresponding research on 
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dust waves have been also extended to the strongly coupled [6-9] as well the magnetized 
dusty plasmas [10-13]. In this letter we study the instability of longitudinal dust wave 
in the strongly coupled dusty plasmas. The stationary external magnetic field is applied. 
Some important factors such as the high-correlation on dusts, the wave-number dependence 
of the elastic relaxation time and effect of dust charging relaxation are included. It is 
shown that there are some frequency regime for the occurence of longitudinal dust wave 
instability because of magnetized effect and the charging collisions between the dusts and 
the plasma particles. The regime of wavenumber that longitudinal wave is allowed depend 
also on the competition between shear viscosity and applied magnetic field strength. The 
behaviors of dust wave instability is also modified significantly in magnetized dusty plasmas 
by magnetized effect. 

The applied external magnetic field is assumed along z-direction as B =Bz. Because we 
consider the longitudinal waves, we just need one of the Maxwell equations, i.e. Poisson 
equation, as 

V 2 </? = -A-ne(5rii - 5n e - Z d0 5n d - n d0 5Z d ), (1) 

where ip is the electrostactic potential, n a o and 5n a are the steady and perturbation density 
of species a = e,i,d (electrons, ions and dusts respectively), Z d0 and 5Z d are the steady and 
perturbation charge numbers of dust negative charge as q d = —eZ d0 — eSZ d . Furthermore, 
the overall charge neutrality condition holds in the steady state, n i0 = n e0 + Z do n dQ . For 
longitudinal modes, one has to assumes that the wave propagates nearly perpendicular 
to the magnetic field. For simplicity we shall consider wave along x. The perturbations 
associated with the problem are assumed as in the form of ~ exp[i(kx — ut)}. Then the 
perturbed electrons (ions) number density in magnetized plasmas can get by ordinary fluid 
equation. On the other hand in terms of a generalized viscoelastic hydrodynamic model [14] 
the perturbed density for the dusts in magnetized plasmas are also obtained. In addition 
to the dust charging equation by orbit-motion-limited (OML) theory the problem would 
be closed by Eq.(l). In Ref. [13] we have obtained the general dispersion relation for low 
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frequency longitudinal dust waves in strongly coupled magnetized dusty plasmas (please 
refer to Eq.(ll) of Ref. [13]) 



where where u; peiPi = ^ Airn e0 , i0 e 2 /m eji , n e0:i0 , m e ,j, V Te ,Ti = \jT e ^/m e ^ and Q eji = eB/m eii c 
are the chracteristic oscillation frequencies, unperturbed number densities, masses, thermal 
velocities and the gyro-frequencies of the electrons and ions, respectively, 

A ~ i + P + l z (l + r + z) [6) 

and 

8=1+ ? T + '\ (4) 

Z(l+T + Z), V 1 

are two parameters which are associated with the dust charging, and 

_ y l fpd /gN 

W w[w + iT7di(A;,a;)] - -f d fi d k 2 V^ d - fl 2 d /[l + rj d2 (k,uj)] 

is the dielectric function of the dust fluid. The other quantities are P = Z do n dQ /n eQ , r = 
Ti/T e is the ratio of ion temperature to the electron temperature, and z = Z d0 e 2 /aT e is 



the normalized dust-dust interaction potential energy, co pd = JAirn d oZ do e 2 /m d , n d o, m d , 



Vrd = \jT d /m d and Q d = Z do eB/m d c are the chracteristic oscillation frequency, unperturbed 
number density, mass, thermal velocity and the gyro- frequency of the dust, especially two 
quantities 

(C + 4q/3)fc 2 

T] dl (k,Uj) = : r, 6 

m d n d0 {l - wr m ) 

and 

m(k,uj) = ; -, (7) 

m d n d0 {l - iujT m ) 

are associated with the bulk viscosity £, the shear viscosity 77 and the viscoelstic relaxation 
time r m for highly correlated dusts. Note that ^ d is the adiabatic index and /i d — 1 + 



u(T)/3 + (T /9)dru(T) is the compressibility and, w(T), the normalized correlation energy, 
or the excess internal energy, of the system [14]. For example when Coulomb coupling 
parameter T = (Z^e) 2 /a d Td, (—Z d e is the dust charge, ad the interdust distance and T d the 
dust temperature) becomes high enough, the negative dispersion occurs for dust waves [9].. 
For weakly coupled dusts (r < 1), or under certain conditions also for strongly coupled dusts 
we have u(T) ~ — v / 3r 3 / 2 /2 [7,8,14]. In practice u(T) is usually obtained by fitting results 
from experiments and/or Monte Carlo or Molecular Dynamics simulations. Accordingly, 
fitting the internal-energy results from a liquid simulation one obtains u(Y) — 0.90r fori < 
T < 200. Other quantities such as the transport coefficients C(r), f/(r), and r m (T) are usually 
too complex to be expressed analytically, so that often simplifying models are employed [14]. 
Introducing the dimensionless quantities u = u/u p d, &d = ^d/^ P d, k = k\ p (where \ p = 



(A, 2 +A e 2 ) x l 2 the plasma Debye scale length, \ = ^Ti/47m i0 e 2 is the ion Debye length, A e = 
^T e /4im e0 e 2 is the electron Debye length), r m = r m u pd , r}{ = (( + 4?y/3) / m d n d QUJ pd a d , rf 2 = 
rj I 'rndndQU) p da d) and \ e ,i,d, P = K,i,d, P /^ P , we can express the general dispersion relation of dust 
waves Eq. (2) as 



2 ,2,2 , iujk 2 K 2 <ql Q 2 

w - <y d /J>dk \ d + 



1 — iuJT m 1 + ik 2 K 2 r]2/ (1 — iuJT m )uj 
1 



1 + £(£;)' 

where k = ad/X p is the ratio of interdust distance to plasma Debye length and 



(8) 



with 



n 2 n 2 

U p i U pe 

For weak magnetized plasmas or unmagnetized plasmas we have C(k) ~ A/k 2 and Qd 
1 so that Eq.(8) will be reduced to the former extensively studied cases [7,9]. However for 
relatively strong magnetized plasmas we have 



here we ignore the contribution from electrons since usually A e ^> Aj and .Q e <C u pe . In 
the following we would study two typical cases, i.e., the kinetic regime of ujr m <C 1 and 
hydrodynamic regime of ur m ^> 1. 

In the kinetic regime of ujr m <C 1 from Eq.(8) we have 

Furthermore we discuss two limiting cases. 

(a) For long-wavelength wave, i.e. the wavelength is much longer than the inter-dust 
distance, kn <C 1, we have 

' + " V (" : + 3«) - n * - ^ ' MTU** 2 -vfT3- a (12) 

where = VLi/uj pi is the normalized Alfven velocity by light speed c. From this we can 
always get k = k r + iki that means there exist wave instablity. Now we look when we can 
get the pure exponent-type unstable wave which needs k r = 0. In fact the frequency of 
wave is in the same order of or less much than the dusty plasma frequency uJ v d so that the 
normalized uo < 0(1) or <C 1. On the other hand Vl d < O(uo) or <C oo and rf 2 ~ 0{j]\) or 
< Thus if Vj\ ~ 0(1) or < 1 then the second term of left hand in Eq.(12) can be ignored. 
For this case where the shear effect is not strong we get 

^ - ng - vlj{vl + A) 

A/(Vl + A) 2 + ^(l-0AT)/3V 1 ' 

Usually Qj <C u pi and 



,>,^^ + ^ (14) 
that is a hybrid waves of dust and ions in the value of less than the dust plasma frequency. 



For r < 1, z ~ 0(1), in general A « 1 to 2 then it leads to that when k > k c = y 15/2.4. 
and T > r c = 2.5k 2 / (k 2 — ft 2 ) then /c 2 < that indicates the existence of pure imaginary 



part of wavenumber. For example for A = 2, k c rs 1.9365, then for k = 2,3,4,5 we have 
correspondingly T c ps 40,4.3,3.3 and 3 respectively. 

Now let us disscuses the condition of ?y* < 0(1) and ujr m <C 1 for the limitation of Y. In 
general it is difficult to obtain exact expressions for r m , rji and r^ 2 etc. However, approximate 
models can be constructed from the results of theories and numerical simulations. For 
example, for the range 10 < T < 200, one can obtain from Table 5 and Fig. 33 of Ref. [14] 
the dependence r]{ ps 0.02a/T from a one component plasma theory. On the other hand 
in terms of r]l/r m = K~ 2 \ 2 d (1 — ^d^d + 4w(r)/15) we have r m ps 3y/T/8. Obviously in this 
range r/l < 0(1) and wr m <C 1 are both hold for ujh < 0J € 1. 

Now we have concluded that in this case, long-wavelength wave in the kinetic regime 
of ujT m <C 1, the spatial instability occurs for moderate and high strong-coupling, weak 
viscosity and moderate viscoelasticity of dust, hybrid order of low- frequency of wave and 
moderate inter-dust distance. 

(b) For short-wavelength wave, i.e. the wavelength is in the same order of or less than 
the inter-dust distance, kn < 0(1), and ^ ^ O(r^) 1, we have 

u 2 -l d fi d k 2 \ 2 d -l = (15) 

and therefore 

e = ^(i-o.4r)/3r- (16) 

That means the wave instablity can occur only if T < 2.5 for u < 1 and 10 > T > 2.5 for 
uj > 1 and 0(1). In later case T < 10 is need for meet the conditions of weak shear viscosity 
rjl 0.02\/T < 1 as well r m ps 3VT/8 0(1). The wave instability occurs independently 
to the dust changing relaxation and the applied magnetic strength. Certainly Vtd/^ P d «C 1 
should not be breaked. 

Now we have concluded that in this case, short-wavelength wave in the kinetic regime of 
wr m 1, the spatial instability occurs for moderate strong-coupling, very weak viscosity, 
weak viscoelasticity of dust, dusty-plasma characteristic oscillation order of frequency of 
wave and has nothing to do with inter-dust distance. 
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In the hydrodynamic regime of ujr m > 1 we have 

^ ~ kV T m ~ r-vZw ~ ^ -T^cWr - (17) 

Similarly we discuss two limiting cases. 

(a) For long-wavelength wave, i.e. the wavelength is much longer than the inter-dust 
distance, k< 1, Eq.(17) would be reduced as 

- 2 - kV i - 4 + &) - ^ - WZTat* -vBa^ < 18 > 

We know that In the long-wavelength limit the k dependence of r m may be modeled by 
[14] rfjT m = K~ 2 \ 2 d (1 - ldf i d + 4«(r)/15). On the other hand (ft 2 ,/^ 2 )^*/^) ~ O^/rJ. 
Therefore one can get the wavenumber of dust wave as 

k 2 „ ^-nl-vl/(vl + A) 
~ a/(vI + a) 2 + k 2 (i-o.3F)/t 1 ; 

which depends on magnetic field and other parameters. Similarly in the case f2j <C u pi and 



u>u H ^^ 2 + ^^ (20) 

we have, therefore, when k > k c « y/lO/3A and r > T c « 3.3k 2 /(k 2 - k 2 ), then k 2 < 0. 
For example for A = 2, k c 1.3, then for /c = 1.4,2,3,4,5 we have correspondingly 
T c « 22.0, 5.7, 4.0, 3.7 and 3.5 respectively. 

It is worthy to note what does mean ooT rn 3> 1. In general u <C 1 or ~ 0(1) so that 
r m ^> 1 is at least required. That means the highly strong-coupling is needed. 

Now we have concluded that in this case, long-wavelength wave in the hydrodynamic 
regime of ur m ^> 1, the spatial instability occurs for highly strong-coupling, highly strong 
viscoelasticity of dust, hybrid order or/and dust-oscillation order frequency of wave and 
moderate inter-dust distance. 

(b) For short-wavelength wave, i.e. the wavelength is in the same order of or less than 



the inter-dust distance, kn < 0(1), and, when Vt d <C yV2/ T m for moderate or/and strong 
shear dusty plasmas, we have 



u 2 -k 2 K 2 ^- ld ^ d k 2 \ 2 d -l = (21) 

or equivalently 

,2 = ^ 2 ~ 1 (no) 

/t 2 (i-o.24r)/3r 1 ; 

that shows the wavenumber k 2 < for T < 4.16 when uo < 1 and otherwise T > 4.16 when 
uj > 1. Obviously the former case is not realized for the requirement of r m ^> 1. Certainly 
the later case holds for very highly strong coupling dusty system. On the other hand, when 



^ yV2/ T m for very weak shear dusty plasma, we have not get the existence of pure 
unstable wave because of its contradictive to the requirement of ur m ^> 1. 

Now we have concluded that in this case, short-wavelength wave in the hydrodynamic 
regime of ur m ^> 1, the spatial instability occurs for highly strong-coupling, moderate 
or/and strong viscosity, highly strong viscoelasticity of dust, dusty-plasma oscillation order 
frequency of wave and has nothing to do with inter-dust distance. 

Finally let us discusse briefly that the condition for the occurence of wave instability in 
the regime of ur m ~ 0(1). In this case the occurence condition of wave instability is almost 
same as in the kinetic regime of UT m <C 1 except r m ~ 0(1) add limitation to T. In fact in 
the middle regime of ur m ~ 0(1) we have 



^ - k ^ 2 - V - i ^ V ^ 2wi -wo - - = 0- (23) 

r m 2 1 - k 2 K, 2 {rjyr m uj 2 )(l - 1)/2 1 + C{k) 

For very small kn <C 1, we have 

,,2 1 _ *Z~2,2 ( Vl | V2 1 r>2 „, ,, J.2 \2 A ,2 n (1A\ 

since k € 1 and r]{/T m 0.05 so that the second term is ignored in Eq.(24) and it is 
recover to the Eq.(13). For large kn ~ 0(1), we have 

^ - ¥ kv i + f 1 + - ^ - 1 - °- < 25 > 



When rjl/r rn <C 1 and <C \/V2/ T m, it leads to 



co 2 - ldl i d k*\\ -1 = (26) 

which is same as Eq.(15). However ur m ~ 0(1) indicates that the occurence condition of 
wave instability in this case requires that the wave frequency is moderately increased and 
the strong-coupling is moderately decreased, compared to the case of ujr m <C 1. 

In summary, we have investigated the low frequency longitudinal wave instability in 
strongly coupled dusty plasmas in presence of magnetized field. The effect of dust charging 
relaxation is taken into accoun. It is found that the occurence of instability of wave is 
determined significantly by the frequency of perturbation wave, the dust charging relaxation, 
the shear viscosity and viscoelastic relaxation time, the coupling parameter of high-correlated 
dusts as well the strength of magnetic field. 
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